o 
o 

(N 



Q 

en 
(N 



X 



hep-ph/0009102 
DESY 00-132 



The (7* -^ qq) - Reggeon Vertex in Next-to-Leading 

Order QCD 



J. Bartels, S. Gieseke, C.-F. Qiao 

fiir Theoretische Physik, Universit; 
Luruper Chaussee 149, 22761 Hamburg, Germany 



CD ■ II. Institut fiir Theoretische Physik, Universitat Hamburg 



> 

(N 

O ■ Abstract 



ON 

^^ I As a first step towards the computation of the NLO corrections to the photon impact 

^P ' factor in the 7*7* — > 7*7* scattering process, we calculate the one loop corrections to 

the coupling of the reggeized gluon to the 7* -^ qq vertex. We hst the results for 

i-^ I the Feynman diagrams which contribute: all loop integrations are carried out, and the 

I ' results are presented in the helicity basis of photon, quark, and antiquark. 

Oh! 



PACS number (s): 11.55.Jy, 12.38.Bx, 13.60.-r 



1 Introduction 



The experimental test of the BFKL Pomeron |I|] is generally considered to be an important 
task in strong interaction physics. Recently much interest has been given to the total cross 
section of the scattering of two highly virtual photons a^^^ 0, ^. This process describes 
the scattering of two small-size projectiles, and its high energy behavior (at not too large 
energies) is expected to be described by the BFKL Pomeron. Therefore, a measurement of 
the reaction e"'"e~ —>■ e^e~ + X hy tagging the outgoing leptons at LEP or at a future Linear 
Collider provides an excellent test of this very important QCD prediction. 

So far leading order calculations of the BFKL Pomeron have been compared to LEP data 
(both OPAL and L3) |Q, ^, |^ . In both experiments the data lie above the one gluon exchange 
curve (commonly called Born approximation), but below the BFKL prediction. Since the 
next-to-leading order (NLO) corrections to the BFKL kernel have been calculated [|^, ^, 
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Figure 1: Contributions to the photon impact factor. 

it is known that the higher corrections will lower the theoretical predictions of the cross 
section. However, a consistent comparison with the NLO BFKL calculations has not yet 
been possible: there remains the task of calculating also the next-to-leading order corrections 
of the coupling of the BFKL Pomeron to the external photons, the so-called photon impact 
factor. 

The photon impact factor is obtained from the energy discontinuity of the amplitude 7* + 
reggeon — * 'y* + reggeon (Fig. 1). In leading order as this discontinuity is simply the square of 
the scattering amplitude 7* + reggeon -^ qq in the tree approximation, and the reggeon, i.e. 
the reggeized gluon, can be identified with the elementary t-channel gluon (with a particular 
helicity). In the next-to-leading order new contributions have to be calculated. For the qq 
intermediate state we need the NLO corrections to the 7* -|- reggeon -^ qq amplitude either 
on the Ihs or on the rhs of the discontinuity line, and the qqg intermediate state requires 
with leading order amplitudes 7* -|- reggeon — > qqg on both sides of the energy discontinuity 
line. The task of calculating the NLO corrections to the photon impact factor therefore 
can therefore be organized in three steps, (i) the calculation of the NLO corrections to the 
7* -|- reggeon -^ qq vertex, (ii) the vertex 7* -^ qqg in leading order, and (iii) the integration 
over the phase space of the intermediate states. In this paper we report on results of the 
first step, the NLO corrections to the '-/* + reggeon -^ qq vertex. The vertex is obtained from 
the high energy limit of the scattering process '~f* + q ^ qq + q. 



2 Technical preliminaries 



The kinematics is illustrated in Fig. 2. As usual, q and p denote the four momenta of the 
photon and the incoming quark, resp., and eL,t the polarization vectors of the photon. We use 
s to denote the energy of the 7*g scattering process process, and we introduce the invariants 



Q' 



-g^ ta = F, tb = {q- k 



M' 



{q + r) , t = r , and x = Q /2p ■ q for the 



Bjorken scaling variable. For simplicity, in the calculation of this paper we treat the quarks 
as massless. The momenta k and r can be written in the Sudakov decomposition form, i.e., 

k = aq' + pp + k^, (1) 




Figure 2: Kinematics of the process 7* + g — * gg + g. 



t , ta + th 

r = -q p + r± , 

s s 



where q' = q + xp and 



f3s 



1 — a 



-Q' 



(2) 



(3) 



The Feynman diagrams which contribute to our NLO-calculation are hsted in Fig. 3. In 
addition to the graphs shown (all diagrams, except for Fig. 3.14), we have to add those 
diagrams where the t-channel gluon couples to the outgoing antiquark rather than the out- 
going quark. It is easy to see that, for the color octet t-channel configuration, the sum of all 
diagrams has to be antisymmetric if we interchange quark and antiquark: k ^ q — k — r, 
A — s> A'. In particular, the 'box' graph shown in Fig. 3.14 has to be antisymmetric by itself. 
We will use the Feynman gauge throughout the calculation, and for the t-channel gluons we 
decompose the metric tensor according to 



9t^u = -{Pf^ql+Puq'J+g^ 



pbV 



(4) 



In our calculation we retain only the first term, since the remaining ones are suppressed by 
powers of the energy. We use the helicity formalism, and our results will be expressed in 
terms of the following matrix elements: 



Tja 
rlrp 

fja 
rlrp 

h: 

Tja 



u{k + r,\) ^ }i ^X%{q-k,X') , 
uik + r,\) ^{4- }i- f) ^X^v^q 
u{k + r,\) ^X'v^q-k^X') , 
u{k + r,\) ^X%{q-k,X') , 
u{k + r,X) ^X%{q-k,X') , 



- k, A') 



(5) 
(6) 
(7) 
(8) 
(9) 




(4) 






(12) 





Figure 3: Feynman diagrams for the process '-/* + q ^ qq + q. 



where A" fa 



1, ■ • ■ ,8) are the color matrices, and A and A' denote the hehcities of the 
outgoing quark pair. When interchanging quark and antiquark hues it will be convenient to 
use the following identity 



H^ + H^ = sH^ - 2e ■ pH^ - 2e ■ tH; . 



(10) 



We organize our calculations in the following order. We first consider those diagrams 
(Fig. 3.1-9) which can be viewed as "elastic scattering of two quarks", with the upper 'incom- 
ing' quark carrying the mass ta- Correspondingly, the diagrams not shown in Fig. 3 define 
quark-quark scattering with the upper 'incoming' antiquark having the mass tj,. Together 
with the correction at the photon vertex (Fig. 3.10) and the quark self-energy (Fig. 3.11), 
these diagrams are the ones which have to be made ultraviolet finite by renormalization. In 
the final part we turn to the calculation of the box diagrams in Figs. 3.12 and 3.14, and the 
pentagon graph shown in Fig. 3.13. 



We are interested in the high energy limit, where 



;ii) 



and we do not impose any restriction on the remaining invariants. The Regge ansatz for the 



scattering amplitude 7*g ^ (99) Q' takes the form: 



rp -r^a _ 



3f^B 



"pO, 



(12) 



Here, 1 + a; is the gluon trajectory. Expanding all terms in powers on the strong coupling g, 
we have 



UJ 






qq 






qq 



qq 



(13) 
(14) 
(15) 



After substituting corresponding elements in Eq.(0) with above expansions, one obtains the 
following structure of the amplitude 



T = /t(°) + (7'T« 



(16) 



with 



and 



7^(0) ^ p(0),a ^ p(0),a 

I'^qq f 11 



7^(1) ^ p(l).« ^ p(0),a , p{0),a ^p(l),a 

I'-^qq f 11 i*^qq +11 



^ 7 ^99 ^ 



1 -s , — s 
m h m 

-t -t 



p(0),a_ 



gq 



(17) 



Here, the e/ represents the charge value of the interacting quark, and the lowest-order 
expressions on the rhs of ([l^) are 



^T = -«(P - r, A,0 /AMp, A,) 



qq 

,(1) 



co^^^it) 



p(0),a 
^ ■y*^qq 



with D = 4-2 e and 

r(l + e)r2(l-e) 
cr = 



(47r)2-^ , I ^^ ' 

(Tja fja 
rirp rirp 
i 
Sta Stb 



r(l-2e) 



7r 



l-^Ee+-{YE-Y)' +^(' 



(19) 
(20) 

(21) 
(22) 



Here Ag and Ag' are the helicities of the incoming and outgoing quark, resp., and A" are the 
generators of the colour group. In this paper we will present the results of Figs. 3.1 - 3.14 
which can be cast into the form of the rhs of (0). AH pieces except for rL_J^^ are known 



from earlier calculations. In particular, the Born result for Vl^J^^gg is well known in the 
context of the photon wave function formalism, and it is available for definite helicity states 
0. The higher order corrections to the quark-quark- reggeon vertex, g^Tqq{t), have been 
calculated in |T^ . For vanishing quark masses these corrections are (following the notations 



of PI): 



r«''^(t) = A" (r«(gg - state) + T^^J{gg - state)) 6, 



AqA^/, 



with 



rjjn?^ - state) 
r«(^(7- state) 



-ty 



(47r)2- 

NJ- 



2 ed 



cr 



ty 



(4vr) 



-2e) 
Icr 



3-2e 



+ 



■^ + — 
e2 3 e 



13 

18 



71' 



2Nr 



.(l),a 



2 3 

- + - + 
e^ e 



(23) 

(24) 
(25) 



Inserting these results on the rhs of (|l8l) we can easily obtain F!^*^^^, which is the goal of 
this paper. 



For future purposes it will be important to note that the vertex FL'^1^ = in (|T2D is expected 



(l),a 



rT| to have a rather complex structure. For example, in the limit of large diffractive masses 



it will be more convenient to write (|T2D as a sum of two different expressions: the first one 
depending on M^ and s, the second one on as and s (or (1 — a)s). This decomposition 
exhibits the reggeization of both the gluon and the quark. A detailed discussion of the 
large-M^ limit will be presented in a seperate paper |12]. 



3 Analytic Results 

We write the NLO amplitude T^^^ for the process 'j* + q 
Feynman diagrams: 



13 



T(i) = 5^ (A, + A,) + A,,. 



gg + g as a sum of the different 



(26) 



j=i 



The subscripts i refer to the numbering in Fig. 3, and the amplitudes Ai correspond to the 
diagrams which are not shown in Fig. 3: they are obtained by interchanging the couplings 
of t-channel gluons between quark and antiquark lines. Formally, we substitute 



Tja 



q — k — r, a 



:i~a] 



) t-a 



th, 



-e- (k + r), X ^ A', 



(27) 



Under these replacements the matrix elements ([7p^(0) remain unchanged. In addition, for 
all diagrams except for Fig. 3.14 we have to include an overall minus sign when interchanging 
quark and antiquark. A14 is already antisymmetric by itself with respect to the interchange 
of quark and antiquark. 



6 



3.1 Calculational methods 

Before presenting explicit analytic expressions for all diagrams we briefly outline the methods 
we have used to obtain the results. 

Starting from the standard Feynman rules of QCD we project the color in the t-channel into 
the antisymmetric octet and proceed by introducing Feynman parameters Xi to combine the 
denominators for the integration of the loop momentum. Only for the simplest diagrams 
these steps are easily performed by hand. Particularly for the diagrams Fig. 3.12-14 this 
step becomes too tedious. Therefore, we have used the computer algebra system Mathe- 



matica with the package FeynCalc |jT5[, in order to reduce the numerators to expressions 
which contain the helicity matrix elements (^ to (P), monomials XiXj ■ ■ ■ in the Feynman 
parameters and, of course, the kinematical invariants. Those diagrams, where the loop in- 
tegral itself does not depend on the large scale s, the box diagrams Fig. 3.12 and Fig. 3.14, 
have to be calculated exactly. The only high energy approximation results from the gluon 
nonsense helicity (i.e. the first term in (0). Diagram Fig. 3.13, on the other hand, has first 
been calculated exactly, and then the high energy limit (0) has been taken. 

To consider the loop integrations, we had to deal with integrals over the loop momentum 
itself, which could be easily performed by the usual shift, and with the remaining integrals 
of the Feynman parameters. Integrals of the latter type are either known, or they could be 
obtained from known ones with the help of recurrence relations to 0{e^) in dimensional regu- 
larization. The technical background for these methods is given in |jl3[ and has been applied 
to almost all cases we are interested in by ]14|. We have written a Mathematica package to 



easily call the results of |T^ . In most cases the integrals are recursively expressed in terms of 
special functions, representing a particular combination of logarithms and dilogarithms for 
a given ra-point function. In the case of diagram Fig. 3.14 we have calculated explicit results 
for integrals with two and three Feynman parameters in the numerator using the derivative 
method |[T^. Also for this task we have used Mathematica. 

Finally, after carrying out the integrals, i.e. replacing the monomials of Feynman parameters 
in our amplitudes with the explicit expressions from our Mathematica package, we have used 
again FeynCalc to take the high energy limit in the diagrams Figs. 3.2, 3, and 13, and to 
carry out some simplifying algebra. A few final simplifications had to be done by hand. The 
expressions that we have obtained in this way will be listed in the remainder of this section. 



3.2 Quark-quark Scattering 

It is suggestive to view the diagrams Figs. 3.1-9 as a quark-quark scattering processes with 
one of the incoming quarks being off-shell (with virtuality t^). We focus on those diagrams in 
Fig. §, where the t-channel gluon(s) couples to the quark. Those diagrams where the gluon(s) 
couples to the antiquark are easily obtained by performing the substitutions described after 



Diagram Fig. 3.1 is one of many three point functions with two massive external legs. Loop 
integrals of this type are well-known, and we have calculated them both by hand and by 
using computer algebra as outlined above. 



Ai 



Nc 1 
2 (47r)2-^ 



^(0) 



2cr tg 
e2 t-U 



cr, , , 2t t 



{-ieef)a 2s 



:fja._-r{0),a 

t-U ' t "" 



X 






* -Ini 



t- Cq t(J 



Here and in the following, we use 



^^°^ = --,^|rS)- 



as a shorthand notation. 



(28) 



(29) 



In the diagrams Fig. 3.2 and Fig. 3.3 we have to deal with four-point integrals, where one 
of the external legs is off-mass shell. Here we make use of the shorthand notations listed in 
the appendix. The result for the sum of the two diagrams is: 



A 



2+3 



- A^iV, 



1 rcp 



t-K 



(47r)2-^ I e2 
ik-t)-' - (- 



{asY' + {-asY' + (-t)-^ - 2[-taY' 
taY") 1 + LdJ'"(as, t, ta) + LdJ™(-as, t, t„) 



with 



Ld^'"(«s,t,4) 



T • /. ^ \ T • /1 (^^\ , t ^ as n^ 
Li2(l - -) + Li2(l - — ) + In-ln- — . 



U 



tr, 



^a ^a 
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Here, Li2 is the standard dilogarithm function defined as 

f^ ln(l - xt) 



Li2(a;) 



dt . 



In order to exhibit the energy dependence we rewrite this expression: 

iV. r cr(-t)- 



A 



2+3 



^(0) 



(47r; 



2-e 



, as , —as 
In h In 

-t -t 



t 



+ 2Li,|l-- + In^ 



J- \ 2 2 



^^H-ty 



2c 



r ^2 



+ 



cr 



i-ty 



2(-U 



t-ta 



{{-tr - {-taY') 



(30) 



(31) 



(32) 



From these two diagrams (and from their respective partners A2+3) we already get the 



complete Ins dependence of our result, in agreement with the rhs of (18). In detail, the first 



line of ( |32[ ) contains, apart from the In a dependence, the leading order trajectory function 
u)^^^{t) (cf.(|20|)). Furthermore, if we take the limit t^ — > in the last two lines, we are left 
with only the last term of the second line. One half of this term will go to the lower vertex 
(|25|) (the same result could also have been deduced from fl^l)- Therefore, the contribution 
from diagrams 2 and 3 to r!^*^^^ is just 



zee/ 



E^ 



Nr 



Sta [Att] 

-ty 



2-e 



2cr(-t)- 



2c 



— 71 



+ 



e 

cr 



In a + 2 Li2 



i-i 

tn 



+ 



In — H 

ta 3 



;-t)--2(-t,)-^- 



t-u 



ii-t)-' - i-ta)-') 



(33) 



We emphasize, once more, that the two diagrams A2+3 + A2+3 provide the complete In s 
dependence on the rhs of (p!2D. At first sight one might expect that also the pentagon 
diagrams Ai^ + Ai^ might contribute to the energy dependence. Later we will show that, in 
fact, this is not the case. 



The calculation of Fig. 3.4 is very similar to the calculation of Fig. 3.1. We find: 



A, 



1 



2iV, (47r)2- 



^(0) 



2cr 



t 



e' t-ta 



ii-ta) 



-tri 



_4 + ^i±^lnl 



,2s. 



-zee/aif^-rW" 



2crt((-ta) 



{-ty 



-ty 



{t-t,y 



2cr(-t)-' 3t + 24, t 5 
h -. r^ln 

€ I Z(i [Z Zq^j Zq^ Z Za- 



(34) 



Diagrams Fig. 3.5 and Fig. 3.6 are needed to obtain the NLO corrections to the lower vertex 
Fqg . In agreement with |10| we find: 



A^ = -A(o) 



NcCr 
(47r)2-^ 



-ty 



2e 



and 



^ (47r)2-^2iV/ ' 



+ - + 



(35) 



(36) 



We note that they can be obtained from Eqs. (^) and (Q) by choosing e < and taking 
the limit t„ — ;> 0. 



The calculation of diagrams Figs. 3.7-9 is straightforward. They contribute to both the 
upper and the lower vertex: 



A 



7+8 



^(0) 



NcCr 



-ty 



(47r)2-^2e(l-2e) VS - 2e 



1 



(37) 



Similarly, the gg-contribution to the gluon- self-energy in Fig. 3.9 leads to 

i-t)-' f 1 



Ap 



^ (47r)2-^e(l-2e) V 3 - 2e 



-1 



(3^ 



These diagrams contribute with equal weight to both the upper and to the lower vertex. 
Therefore, in order to complete the upper reggeon- quark- quark vertex with one off-shell 
quark, we simply add 1/2 of the sum of (0) and ( pS]) to the contributions (pSf), (0), 
and (|35|) . Similarly, the lower reggeon-quark-quark vertex with all quarks being on-shell is 
obtained by adding 1/2 of the sum of (pTf ) and ( |38D to (|35|), ( |36D and the contribution of the 
box diagrams Figs. 3.2 and 3.3: 



.^(0) 



N.. 



-ty 



(4vr; 



2-e 



2cr 



TT 



(39) 



To summarize, so far we have analyzed the diagrams Figs. 3.1-9. From Figs. 3.5, 3.6, 1/2 
of Figs. 3.7-9, and from a piece of Fig. 3.2 and 3.3. we have reproduced the lower reggeon- 
quark-quark vertex of [0. From Figs. 3.1, 3.4, 1/2 of Figs. 3.7-9, and from the major part 
of Figs. 3.2 and 3.3 we have computed a new reggeon-quark-quark vertex with one quark 
having the mass ta. Moreover, from Figs. 3.2 and 3.3. we have extracted the Ins terms of 
the rhs of (|TB|). The remaining diagrams Figs. 3.10-14 provide contributions to the upper 
vertex Tl^J'Z^q only. 



3.3 Vertex correction and quark self-energy 



In this subsection we present the results of the vertex correction Fig. 3.10: 



A 



Cf 



10 



(4vr) 



2-e 



^(0) 



2cr g' 
e2 Q-^ + ta 



{{Q 



2\-e 



cr 



i-taD+-i-ta 



(40) 



+ 4 + 



3Q' 



In- 



-U 



Q2 + ta Q 



t '^'^ 



with 



^f££H„,-.,io^i.^^Ji„^, 



r = —lecf 



4cr Q^ 

e^ Q' + ta 

t, 



{m- - i-taD 



and for the quark self-energy Fig. 3.11: 

^1 = -A(o) 



Cf 



xr(l-e) 



-t 



(47r)2-^' "' e(l-2e)' 
Eqs. (|40| ) and ( ^2] ) provide new contributions to the vertex T^^Jf^^^. 



(41) 



(42) 
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3.4 The Box diagram Fig. 3.12 

In this subsection we give the result of the diagram Fig. 3.12, which was calculated entirely 



with the help of the computer algebra (in Jill ^^^^ diagram has been named 'adjacent box'). 
The result will be expressed as follows: 



A 



12 



r(o),- 



„2s f-ieej 



ct 



'^ t \ s J \ 2NJ (47r)2-^ 
Starting with the divergent terms, we have 



2^12 + ^12 



A 



A 



(-2) 



2aH^st{{~tr~{-tar) 



12 






2h:^ 



-ty 



-M 



2\-^ 



+ {Q 



2\-^ 



t-ta 



(43) 



(44) 



Q^+ta 






{{Q')-' - i-tar) 



and 



A 



(-1) 

12 



-Ae-km 



+ ms 



ii-ty 



-ta) 



i-tay 



M2 Q^+ta 



t-tn 



M2 



(45) 



The 0{e^) term for the this adjacent box is expressed in terms of many different functions, 
related to this loop integral. For convenience, the definitions of these functions are listed in 
the appendix. 



A 



(0) 
12 



-2{2H!^ + saiJ;)Lco(-g^ M^, t) 

{3H^M'^ + H^[Qe ■ pM^ + 2s{{8a - 3)e ■ k + (5a - 3)^ ■ r)] 



Lci(M2,-Q2^t) 



- saH^{3ta - 2U - Q') - 2H;{2Qh -r + e- k{M^ - Q^ - 3h))}- ^^ 

+ 2{H^M^ + 2H^{M^e ■ p - s[e ■ k + {1 - a)e ■ r]) + 2sH^[{l - a)t + aQ^ + 4] 

+ {2e ■ kH;{2M^ - 2Q2 _ tj _ (i _ a)sQ^H'', - H^{Q^ + Q 



AP 



- 2Hl[e ■ k{l - 2a)s + e ■ p{Q^ + t„)]}- 
- 2{2sa(3e -k + e- r)Hl + sH''^[a{M'^ - Q^ - 3ta) + (t - 4)] 

+ H;[3e.rt^ + e.k{At-Q]]^^^S^ 
+ 2{Ae- pM^ HI + 3M^H^ - s[a(M=^ - 2 Q^ - t) + t]H'^ 
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+ 4as/J,"Lc^'"(t„, t)-Se-k H;LcI"'{~Q\ Q 

- 2 {H^ M^ + 2H^[e ■ pM^ - s((l - 2a)e ■ k + {1 - a)e ■ r)] 

+ 2H;e-k{2M' + t-Q-asH:{Q' + Q}^^^^^t^MlA 

+ 2{2£ ■ k{ta - m^)h; + (g2 + ta)H^ 

+ 2[{l-2a)se.k + e.p{Q' + Q]Hl] ^'''^^'^';^''^^ 

- 2(25 ■kH'; + saH^)Lcs{-Q^M',t) - 4 h^^^<^^^1^^1izS1A 

+ 2{tbM^H^ + sta[a {Q^ + t)-Q^-2aM^ - U]H^ 

- 2ta[e ■ k {2M^ + Q^ + t)+e-r{t- h)]H^ - 2[{-{e ■ pM^) + {e ■ k + e ■ r)s)ta 

- 2{H^ + sH:-2e-r H;}Ldis{ta, M\ -Q\ t) 

+ {A[as{Qe ■ kit -U) + 2e ■ kta + 2e ■ rQ"^ + ?>e ■ rta) 

-ta{e-pM^ + s{€-k + e- r))]Hl + A\2e ■ rta{t -h) + e ■ k {2t{Q^ + t) 
+ {t- ta){ta - 2U))]H; - 2 s[a{Q^ + 4) (2 Q^ + 54) - 2{Qh + tah)] H^ 

+ 2{ - (g2 + QM^H^ + 2[e-r (Q^ + Q (t - t,) 

+ e ■ k{{Q^ + ta){t -ta) + it -Uf)]H; -2[e ■ pM\Q^ + ta) 

- s{{l - a)e ■ r {Q^ + ta) - e ■ k{Q^ + ta- {I -^ a) M^))]H^ 

+ s [(1 - a)(g t + tah) -a{Q + to) ]H^] j^^ 

^ 2 S ^ (g^ + 2 4) ff," Ld24(ta, M\ -Q\ t) 

+ {24e ■ pM^ i/^ + 18 M'^H^ - 6 s[t - t„ - a (3 g^ + 24 + h)]H'^ 

- 12[. ■ A: (3 1 - 24 - 4) + . ■ r (g^ + Q )]H;] ^^^^^'"'^^^ "^'' '^ 

+ {8£ ■ A; 4(4 - t)/j;^ - A{QH + 44) iJ^ + 8[5 ■ A; s (4 - t) - a e ■ A; s (M^ - 2 1 + 2 4) 

-e-p(g t + ta4)]i/fe} ^^^^5^ 

+ {4 £ ■ A; (t - tj'i/; - 4 s [(£ ■ A; + £ • r) (t - 4) - a (2 £ • A; + £ ■ r ) (g2 + 4)]iJ« 

Ld3i4(t„,M2,-g2,t) 



-2as(g2t + 44)ij;}- 



M24 



+ 4 (1 - a) . . . * i/J Ld,...(i„.M^-Q^i) _ 4 „,(,., + ,.,) ^. Ld3„(i., A/^ -Q^ t) 



''a ''a 
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3.5 The Box diagram Fig. 3.14 



The result of the diagram Fig. 3.14 (in [|14| named 'opposite box'), can, again, be spht into 
divergent and finite pieces: 



A 



14 



r(o),. 



a 2s f-iecf 



N, 



"^ t \ s ) \2 
The term proportional to cr/e^ reads: 



[A-n] 



2-e 



cr .(-2) cr (-1) (0)- 

9 ^14 \ ^14 \ ^ 

e^ e 



H4 



(47) 



A 



(-2) 
14 



2 £ ■ pi/," (t„ - u)_ ^ {h;^ - H^) {ta + h) 



QH + tah 



+ 



His 



2{t-taf{t-hf (QH + tah 



QH + tah 

a - {1 - a)) {QH + tah) 



X (t^ (t„ + h) + ta h (ta + h)-Atta U) 

- {ta - tb) (2(t - taf{t - 4)2 + (g2 t + t^ t,) (^tatb - t')) 



2i7; 



(t-t„) (t-tb) (QH + tatb) 
+ e-k{t -ta) (tb (tb -ta) +t {ta + tfc)) |^ (-t) 



e-{k + r) {t- tb) {ta {ta -tb)+t {ta + tfc)) 



+ 



-2e-pH^M^Q^ {ta-tb 



+ 



HI M^Q^s {ta-tb 



+ 



(Q' + ta) (g2 + tb) (Q2 t + 4 tfe ) (Q2 + ta) {Q^ + tfc) (Q^ t + tatb) 

{H^ - H^) Q2 (t^2 + t^2 _ ^ (^^ ^ ^^) ^ g2 (^^ + t, _ 2 t)) 



{Q' + ta) {Q^+tb) {QH + tatb) 

'"" {^-A; {Q^ + tb) (g' (ta-tfe) 



2mQ^ 



+ 



(g2 + 4)'(g2 + t,)2(Q2i + i^^^) 

+ Q2 (2 t„2 + tatb- tb^ + t{tb- 3ta)) + ta {ta^ + tb^ - t {ta + tb) ) 

+ e-{k + r) (g2 + ta) (Q^ {tb - ta) + Q^ {2tl + tatb - ta^ + t{ta - 3tfe)) 

+ tb [ta^ +tb^-t {ta + tb)) )] 1 (g')~' 

{H^-H^) {Q^ + M^Q^-ta^) 2e-pH^[ta^ + Q^ {t + ta-tb)] 



{Q'+ta) {QH + tatb) {Q^ + ta) {QH + tatb) 


H^s 


ata{QH + tltb) -tl{t-ta)^ 


{t-taf {Q^+ta) {QH + tatb) 


+ g2 (_^3^^2(^^ +tb)+tta{ta 


-2tb)+tl {tb-ta)+ata^ {t + tb)) 
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2h; 



e-k{t-ta) {Q'^ + M^Q^ {Q'-ta) 



it-ta) {Q^+taY (QH + tah) 



+ Q^ {ta {ta-h)+t {2ta + h)-t^)) 



-t„ 



^ im - H^) JQ* + M^Q^- h^) ^^e-pHl [t,2 + Q^ {t + u- ta 



{Q^ + U) {QH + tah 



His 



(Q2 + 4) [QH + tah] 

l-a)tbiQH + tlta) +tl{t-thf 



(t-4) (Q^ + h) (QH + tat, 

+ Q^ {t^ - t' (ta + th) +tthi2ta - 4) + tl (tb - ta) " (1 " a)tb^ {t + t„)) 



+ 



2h; 



it-t,) {Q^ + tbY (QH + tatk) 

+ QUb + Q^ tb^ + tb^) +e-k {Q^ + tb) {tb^ (t + tb) 

+ Q^ {tb {tb-ta)+t {2tb + ta)-t^))]\{-tb)~\ 



e-{k + r){t- tb) [Q^ + M^ g2 (g2 _ ^^^ 



(4J 



The cr/e - term is rather simple, since we kept only the first terms in the expansion of powers 
like (— 1)~^ while the logarithms are combined with logarithms from the finite term, leading 
to significant simplifications: 



A 



(-1) 

14 



s Hl[a{t - tb) -{l-a){t- to)] 4 H;[e . k {Q'' + tb) + e ■ {k + r) {Q'' + tj] 



it - ta) (t - tb) 



(Q^+ta) (Q' + tb) 



(49) 



Finally, the 0(e°) term for this opposite box diagram reads: 

(0) _ 2 s [a{t - tb) - (1 - a)(t - 4)] H^ 8[e-{r + k) {Q^ + ta)+e-k{Q^ + tb)]H; 



A 



'' it-ta){t-tb) {Q^ + ta){Q^+tb) 

2s[ae-{k + r){Q^ + ta) {t - tb) + {I - a)e ■ k {Q^ + tb) {t - ta)]H^ 



+ 



{t - ta) {t - tb) (Q2 + ta) (Q2 + tb) 



Ld^{ta,tb,-Q\t) 
M'^iQH + tatb) 



{H^ - H^) M' 3{Qh + ta tb) + M\ta + tb) 



+ 2H\ 



+ h; 



■pM^ {ta-tb)-3s {QH + tatb){ae-k + {l-a) e ■ {k + r)) 

r {ta - tb) {2Q^ + g'(4(t„ + tb) - 5t) - 2t{ta + 4) + 2(t„ + tb)^ - 3ta tb) 
+ i€-k + e-ik + r)) (4tg^ + Q^t(9(t„ + tb) - lOt) - 2t^{ta + tb) 

+ 2t{tl + tl- tatb) + 3tatb{ta + tb)) 
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H-s 



{ta - h){2Q^ - 7QH + AQ^ita + h) + 2(t -ta- hf - 3tah 



- 3 [a - (1 - a)] (g2 1 + tah) (2g2 + t, + t^) ] | 



3 H^ Q^ s 



M2 (Q^+ta) {Q^+U)l 
V 



2mQ'' 



M^Q-^ + tafiQ^ + hf 



«(Q'+ta)-(l-«)(Q'+4) 

e-k {AM^ + ?,{Q^ + Q){Q^ + hf 



M2 (g2 + taf (g2 + 1,)2 (g2 ^ + ^^ ^^) 



e-k{Q^ + h)\Qa{Q^ + ta){t-tk 



+ 2(1 - a)4M2) + 5 . (A; + r) (g2 + 4)2 ^g (^ _ ^^ ^g2 ^ ^^^ (^ _ ^j 
+ 2at5M2)j| log(g2) 

£-(A; + r)(M2(t;,-t„) 



(ifj^ — /Jj^) {2t — tg — th) 



lip 



+ 



2 (t - tj (t - tfe) M2 (t - tj (t - h) 

6t(t - t,)) ~e-k {M^th - ta) + 6t(t - t„)) 

t (ta-tfe) {iAt-2ta-2h)Q^-t^ 



H-s 



2M^t-taf{t-hf 



+ 



+ 3t(t, + t,) - 2t2 - 2tl - tah) 

+ (!-«)(- At\ta + tb) + 4t6(4 + hf + t2(7t2 + Qtah + 7t^) 

- t(3t^ + htlu + 54tfe + 3t^) 

+ g'(-6t3 + 7t2(4 + h) + t„t6(ta + h) - tiStl + 44tfe + 3tl))) 
-a{- At%ta + h) + 4tfc(4 + tfc)2 + t\7tl + 64tfe + 7tl) 

- t(3t3 + St^t, + 5tatl + 3tl) 

+ Q^i-Qt^ + 7t\ta + tft) + tahita + 4) - t(3t^ + 4^;, + 3t2))) 

pM^ {t- ta) {t - tb) {ta - tb) (g' t + ta tfe) 

+ e-kst (t-ta) {2M^tb{l-a){t-ta) + 6a{Q^ + tb){t-tb)^) 

+ €-{k + r)st{t- tb) {2MHaa{t-tb) + 6{l- a) (Q^ + ta) {t - 4)^) | log(-t) 

M^ (g2 + ta) {ta-t + Aat) 



Hk 

M^{t-taf{t-tbf {QH + tatb)^ 



H^-H^ 



+ 



H-s 



2{ta-t) ' 2M^t-taf {Q^+ta) 

- Qtait - ta){ta + ag^ - (1 - a)t) 

2 



h; 



M^{t-ta){Q^+taf^ ^^ ' 

+ 6t„£ ■ r (g2 + taY + 6tae-k {Q^+ ta) (g' + t) + 4£ ■ kM\Q^ - ta){ta " t) 
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+ 



Ht 



+ 



+ se-r {2aMH + 6{l -a){t-taf) ta{Q^ + ta Y + s e ■ k [2 M\[at{Q^ + ta) 

- Q\t - taf - 2aQ\t - ta){t + 2ta)] + 6(t - ta)[ta{t " QiQ^ + t^f 



-aM\QH + tl)]) log(-t„) 



Ija Tja 



+ 



H-s 



2{h-t) 2M^t-t,y {Q-^ + U)^ 
- Qhit - h){at ~ tb - {I - a)Q^) + 

2 



-QUe-r [Q^ + uY + QUe-{k + r) {Q^ + U) {Q^ + t) 
+ A€-{k + r)M\Q^-h)ih-t) 



M^ {Q^ + h){t-U-4{l-a)t) 



P 



Hi 



M^ {t - hf {Q^ + Uf {QH + tah)^ 



-e-pM^ (t-h) {Q^ + h) {Q^t + tah) 

+ se-r {-2{l-a)MH-6a{t-tbf) h {Q^ + h Y 

+ se-{k + r){2 M\[{1 - a)t{Q^ + 4)' - Q^t - hf - 2(1 - a)Q'(t - h){t + 2h)] 

+ 6(t - h)[h(t - h){Q^ + hf - (1 - a)M\QH + t^)])l 1 log(-t,). 



(50) 



We have cast all our results for An into a form which manifestly exhibits the antisymmetry 
under the exchange of q and q. 



3.6 The Pentagon diagram Fig. 3.13 

The result for the pentagon diagram Fig. 3.13, is much simpler than one might expect. The 
reason for this is the fact that, in contrast to the box diagrams Fig. 3.12 and 3.14, the 
integrals depend upon the large scale s, and in the high energy limit they can enormously 
be simplified. As before, we seperate divergent and finite pieces: 



A 



13 



p(0),a^ 



-iecf t 



Ct 



2 ) ■(47r)2-^ 



1 



2^13 



(-2) 



A 



(51) 
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The divergent contribution has the form: 



A 



(-2) 



13 
+ 
+ 



Tja 

llrp 



rja 

llrp 



ta h 



(( 



a s] 



tn 



Tja 

llrp 

Tja 

llrp 



Tja 

llrp 



t{QH + tatb) 

tah 
t{QH + tah) [ta 

_ tatb 

' t{QH + tatb) 



h 

fja 

llrp 



+ 



2H^ 



2H; 



tta QH + tah 



+ 



2H'' 

p 



H!^ H^ 



h 



2H^ 



th Q^t + tah 
+ 



e-kQ^ e-{k + r)ta 
e-ktb e-{k + r)Q'^' 



t-h 



Q^ + h 



2H- 



QH + tah 



e-kth e-{k + r)ta 



t-h 



t-tn 



i-ty 



+ 



Q' 



QH + tatk 



m 



Tja 

llrp 



2Q^H^ 

QH + tatb 



e ■ k e ■ [k + r) 



Q^ + ta Q^+t, 



Q 



2-^ 



whereas the finite piece reads: 



A 



(0) 



tatb 



Tja Tja 

llrp llrp 

ta h 



Ld^(t„(l-a)s,-Q2,-as) 



- Ld^oitb, -a s, -Q^ (1 -a)s)- Ld^'"((l - a) s, t, tb) + LdJ"'(-a s, t, ta) 



2m 



Tt ja 



tta 



Ldl{ta,il-a)s,-Q\-as) 



2H^ 

TTb 



Ldo(tb, -a s, -Q^, (1 - a) s) 



+ 



tat 



a l^b 



Tja 

llrp 



Tja 

llrp 



t{QH + tatb, 
M2 {QH + tatb) 



ta tb 



e-k{Q^ + tb)+e-{k + r){Q^ + ta) ^ Ld°P(4,tfc,t, -Q^). (53) 



(52) 



The counterpart Ai^ of the pentagon graph Fig. 3.13 can be obtained in two different ways. 
Either we follow the substitution described after ( p6D and simply interchange quark and 
antiquark q ^^ q. Alternatively, we start from Fig. 3.13 and perform the crossing (s ^ m ~ 
—s) (an additional minus sign comes from the color antisymmetry in the t-channel). In the 
following we present the sum ^13 + ^13. In order to demonstrate the cancellation of the 
In s-dependence we combine the first line of (0) with the first three lines of ( p3D (and their 
counterparts in A13). The results contain a term proportional to Cr/e: 



Tja 

llrp 



Tja 

llrp 



ta tb 



[ln((l-a)s)+ln(-(l-«)s) 
— ln(a s) — ln(— a s) 



(54) 
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and a finite piece: 

tb H^ + ta H^ 



]- ( In((l -a)s)- ln(-(l -a)s)Y -]-{ ln(a s) - ln(-a s))' 



t{QH + tah) 

+ (H-ta) + H-h) - In(-t) - ln(g2)) 

X ( ln((l — a) s) + ln(— (1 — a) s) — ln{a s) — ln(— a s)) 



+ 



Tja 
rirp 

- ( In((l -a)s)- ln(-a s)) ( ln(-(l - a) s) - ln(a s)) 

+ (ln(g2)_2 ln(-t„)) (ln((l-a)s) + ln(-(l-a)s)-ln(as)-ln(-as)) 

^[i(ln((l-a).)-ln(-(l-a).))^-i(ln(a.)-ln(-a.))^ 

- (ln((l-a)s)-ln(-as)) (ln(-(l - a) s) - ln(a s)) 

+ ( ln(Q2) - 2 In(-tb)) ( ln(a s) + ln(-a s) - ln((l - a) s) - ln(-(l - a) s)) 



;(ln(as) -ln(-as))^ - -(ln((l - a) s) - ln(-(l - a) s))' 



+ 2 






Lis 1 + 






Lis 1 






Lis 1 - 



+ Lio 1 - 



+ -ln2 
2 



-h) 



■In^f-t. 



ln(-4)ln- + ln(-t„)ln- 

tb ta 



T 



tt, 






if^ 






(55) 



We can easily observe that the logarithms in s cancel out. Note, however, the nontrivial 
phase structure: as we have discussed after (P3[), such a phase structure is expected, and 
in the double- Regge limit (t, Q'^ <^ M^ <^ s) it will lead to the anticipated decomposition. 
Finally we mention that, when adding the two pentagon graphs, the remaining four lines in 
(|52| ) and the last two lines in (|53| ) are simply multiplied by a factor of 2. 



4 Renormalization 



In this paper the Feynman gauge is adopted in all calculations. In order to regularize 
the singularities we have used the dimensional regularization procedure. At this stage our 
results contain both infrared and ultraviolet singularities: standard renormalization will 
remove the ultraviolet singularities, whereas the infrared singularities will cancel only in 
the complete NLO result for the photon impact factor |TH]. In this section we will perform 
the renormalization, and we will use the modified minimal subtraction (MS) scheme. In 
order to demonstrate the cancellation of the ultraviolet e-poles, we first list the ultraviolet 
divergencies of our diagrams in Figs. 3.1 and 3.4-11 (the box diagrams in Figs. 3.2 and 3.3, 
the pentagon graph Fig. 3.13, and and the box diagrams Figs. 3.12 and 3. 14 are ultraviolet 
finite). Our analysis leads to the following results (deviating from definitions in Eq. (llO) we 



now include the coupling constant 



A'i 



uv 



^(0)^4 ^_^^-.3^^ 



(47r; 



2-e 



euv 2 



(56) 



A^ 



^(0)^4 



/ {-ty 



[An] 



2-e 



euv 2Ar, 



(57) 



A 



^y_ A(V (-t)-^3iV, 



(47r; 



2-e 



^(yy 



(5^ 



aUV 
^6 



^(0)^4 



/ {-ty 



(47r)2-^ e^v- 2iV, 



(59) 



^(0)^4 



AUV 

^7+8+9 - (A^\2~e 



g' i-ty 



(47r) 



euv 3 



-A^. 



-n 



/] 



(60) 



auv 

^10 



A(°)^^ (-to 



(4vr; 



2-e 



Ci/y 



-Cf , 



(61) 



uv 



AU\ 
^11 



(47r)2-^ e^y 



a 



(62) 



The poles in A^ and Ag, coming from the lower vertex correction, as well as the gluon 
self-energy A^ to Ag can be compared to standard textbook results. 

In above formulas g, the strong coupling constant, denotes the unrenormalized one, the bare 
coupling. In order to perform the usual renormalization, we make the replacement: 



9 -^ r. rrr- "^ 9r 



Zi\^ Z-x 



1 - ?iA 



4vr \euv 



7e + log(47r) + 



(63) 



where /9o = (-g-^c — 3^/)) ^'^'^ 7e is the Euler constant. In the MS scheme we have: 



«., 



Zx = 1--^{N, + Cf) 7E + log(47r) 

47r leuv 



1 



Z2 
Z3 



ttc ,5 



^uv 



- 7e + log(47r) 



1 + 



47r'3' 



Nr 



r^/) 



ec/y 



7e + log(47r) 



(64) 
(65) 
(66) 
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for the vertex renormalization, the quark wave function renormahzation, and for the gluon 
wave function renormahzation, resp. Finahy, we add the quark self-energy diagrams for the 
external legs (not shown in Fig. 3): 



^uv _ ^9 ^F 

(47r)2-euv 



A^^ = ^ ^ (p,7) 

quark self —energy /'/i_\2— f , ' V^ ' / 



It is easy to see that in this way all ultraviolet divergencies cancel. In particular, the 
ultraviolet divergences in diagrams Figs. 3. 10 and 3.11 exactly cancel against each other as 
expected, similar to situation in the NLO calculation of 7* — * qq in the e~^e~ annihilation 
process. 



5 Conclusions 

In this paper we have calculated the high energy limit of the process 7*g — ^ qqq in next-to- 
leading order, and from the results we have extracted the NLO corrections to the coupling 
of the reggeized gluon to the vertex 7* -^ qq. This calculation represents the first step in the 
computation of the NLO corrections to the photon impact factor. These NLO corrections 
will allow to perform a complete NLO analysis of the BFKL prediction for the scattering 
process 7*7* — * 7*7* at high energies. 

In this paper we have listed the results of all one loop diagrams. Using dimensional regular- 
ization we have carried out all loop integrations, and our final results are expressed in terms 
of logarithms and dilogarithms. We also show the explicit dependence upon the helicities of 
the photon and the quarks. After renormalization our results are free from ultraviolet diver- 
gencies, but they still contain infrared singularities which will cancel once all NLO pieces of 
the photon impact factor have been calculated and put together. 

We have not yet attempted to combine the contributions from the Feynman diagrams into 
a single compact expression. The results for the individual diagrams are sufficiently compli- 
cated and lengthy, and we found it useful to first list them separately. A closer investigation 
of the sum of all diagrams will be presented in a forthcoming paper. This includes important 
consistency checks as well as investigations of several special kinematic limits. 



Note: shortly before this paper has been completed a short paper (hep-ph/0007119) by 
V. Fadin, D. Ivanov, and M. Kotsky has appeared which reports on a similar calculation 
of the 7* — > qq vertex. However, the results presented in that paper are written in terms 
of one-dimensional, two-dimensional, or even three-dimensional integrals which have to be 
performed. We therefore feel, at this stage, unable to make any comparison between our and 
their results. 
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Appendix 

A Basic functions 



In this section we define some basic functions, which are closely related to scalar integrals 
and are built up from Logarithms and Dilogarithms. The functions are defined in [|1^. For 
convenience we list them here. 

The triangle functions depend on virtualities pf of the external legs i, which we denote by 
p1,P2 and pi = Si2- In the case of the two mass triangle, where the second leg is on-shell, 
P2 = 0, we have only one simple logarithm: 

Lcr(p?,.i2) = ln^. (68) 

pi 



For the three-mass triangle we have the following function 



LCo(Pi,P2'Sl2) 



with the definitions 



ln(a+a-) In ( ^ ) + 2Li2(a+) - 2Li2(a" 

-L (Jj 



(69) 



As = -pt-pt- s'u + W1P2 + 2PiSi2 + 2p^si2, (70) 

± si2 + pi-p?± V-A3 i^,. 

The boxes depend on the virtualities of the external legs (we put p\ = S123) and on the 
invariants S12 = (pi +^2)^, -523 = (p2 +^3)^- In the case of the one-mass box, with only the 
fourth leg having a virtuality S123 7^ we have the function 

LdlTisu, .23, .123) = Li2 (l-^)+uJl-^)+ln^ln^-^. (72) 

V .123/ V S123/ S123 S123 6 
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For the box with the two adjacent legs '"1"' and '"4"' being off-shell, {pi ^ and pi = 
•Si23 7^ 0) we have 

Ldg(.i2, S2s,pI sus) = U, f 1 - — V Li2 (l-^]+]-ln -^ In ^. (73) 

V ■5123/ V 512/2 pfSus Si23 

Finally, in case of the opposite box, where pi ^ and pi = sus 7^ we have 
Ldo^(si2, 523,^25 -^m) =Li2 ( 1 — I +Li2 ( 1 — I - Li2 ( 1 -) 

V ■5123/ V ■5123 / V ^523/ 



-Li2(l-^)+Li2fl-^+ln^ln^. (74) 

V ■5l2'523/ ■S123 ^5123 



The pentagon (in our case it always has one virtual external leg, pf 7^ 0), will always be 
expressed in terms of these box functions, since it is basically calculated in terms of boxes, 
introduced by removing propagators from the pentagon. 



B More Functions 



In the following we give a list of functions, appearing in our calculations. These functions 
appear in the tensor decomposition of the loop integrals and will recursively be expressed in 
terms of the basic functions we introduced in the previous section. 



B.l Functions for the two-mass triangle 

Here we list the functions which are neede in order to express the tensor integrals of the 
two-mass triangle. The invariants are explained above. 



Lc^(p?,si2)=^^^£^^^ (75) 



S12 - PJ 

^2m/J2 _ 

'5'12i — 

■S12 - Pi 

2t ^ 
J2m/J2 "" 

»12i = 

S12 - Pi 



Lcr(p;,.„) = -?H^(iiiH)^ (76) 

Pi 



B.2 Three-mass triangle functions 

In the tensor decomposition of the three-mass triangle the following functions appear (the 
Gram-determinant A3 is given in eq. ([70|) ): 
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Lci(pi, pI, S12) =-r- [ - 2p^ In ^ + {pi + pi - Su) In -^ 
^3 Pi P2 



+ Pi(5i2 +P2- Pi)Lco(Pi, P2> Sl2)J (78) 

LC2(Pi,P2> S12) =77-^ [Pi + P2 - S12 - P2 In ^ + 2^2(^1 - P2 + Sl2)LCi(Pi,P2, S12) 
2A3 P2 

+ P?(-P?+P2 + Si2)Lci(p2,P?,Sl2) -PiP2Lco(p?,P2>Sl2)] (79) 

Lc3(p?,P2, S12) =777^ [ - 2p? + (P2 - S12) \n^ + p%Co{pI,pI, S12) 
2A3 P2 

+ 3p?(-p? + p^ + si2)Lci(p2,p2^ 5,2)] (80) 

Lci5(p?,pLsi2) =-[piLci(p2,P?,Si2)+P2Lci(Pi,P2>Si2)] (81) 

LC25(P?,pL512) =7T-[2p?P2Sl2LCi5(pi,P2,Sl2) " ^ (Pl (^12 + P2 " Pi) l^ ^ 
4ZA3 O P]^ 

+ P2(si2 + P? - P^) In ^ + 2pIpIsi2)] (82) 

P2 

LC35(P?,pLs12) =— -[2piP2Sl2LC25(P?,P2>Sl2) - —{pl{s 12 + pI - p\) \n -^ 

0A3 oU pf 

+ P2(5i2 + pI- pD In ^ + pIpIsu{pI +pI + Su)/2)] (83) 

P2 

B.3 Adjacent Box 

For the box with two adjacent massive external legs we use the abbreviation 

A4 = 2S23 [(Sl23 - Si2)(si2 " pI) - S12S23] ■ (84) 

With this we have the following functions related to higher dimensionial scalar integrals: 

J A r 2 N 2Si2S23 fy ia/ 2 N 

Ldi5(Si2,S23,Pl,Sl23) = T [Ldg (Sl2, S23, Pi, S123) 

A4 

/ 2 Pl^l23 \ / 2 \1 / \ 

+ o (-^123 +Pl- S23 - 2 )Lco(Pi, S23, Si23)j (85) 

^ S12 

2 
T J / 2 \ S12S23 [-523 , S123 , S123 Pi 1 S123 

Ld25(si2, S23, Pi, S123) = — :tx~ i^r^'^ + ^12 In —\n^ 

6L\4 I S23 S12 l p^ 

+ Sl2S23Ldi5(Si2, S23, P?, S123) + (Sl23 + p\ - ^23 ^ )LCis(pi, ^23, S123)] (86) 

S12 

J . / 2 N S12S23 rSi3 S123 , S?2 1 S123 Pi , S123 , S12S23 

Ldss S12, S23,Pi, S123 = — T-T— — In + 777 In T^ In ^- + —-— 

5A4 '- 24 S23 12 S12 24 PI 12 



2»i2„ 
+ Sl2S23Ld25(si2, ^23, P?, S123) + (^123 + p\ - ^23 ^ )LC25(P?, ^23, S123)] (87) 



Sl2 
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In the tensor decomposition of adjacent box integrals the following functions are used: 

Ldi(si2,S23,P?,Si23) = -[Ldi5(Si2,S23,P?,Si23) + Lco(pi, S23, Si 



2 

Ld2i(si2, S23, pI, S123) = [3Ld25(si2, S23,pI, S123) + Lci5(pi, S23, S123 

S12 

-Lci(s23,P?,Sl23) (89) 

Ld22(Sl2, S23,pI, S123) = 2^^^ — [3Ld25(Sl2, S23,pI, ^123) + ^^isipl, ^23, Sl23 



S12S23 
_ 'JlLc'ripl S12) + '-^LC,{S23, pI, 5123) (90) 

S23 S23 

2 

Ld24(si2, S23, pI, S123) = '2— [3Ld25(si2, S23, pI, S123) + ^Cis{pl, -523, S123 



'1, S123) = 2 

S12S23 

2 

- — Lc^(Si2,Si23) + — LC1(S23,P?,S123) (91) 

S23 S23 

/ 12 
Ld3ii(si2,S23,P?,Si23) = " (^^- [5Ld35(si2, S23, P?, S123) + Lc25(Pl, 523, 5123 ) 



S12S23 

i 
S23 ^ ' "^ '"" S23 

,Sl23) = -(^ 

+ — ^LC3(s23,J9i,Si23) + — Lc2(Pi, S23, S123) ) (92) 

-512 S12 / 

2 

Ld3i4(si2, S23,pI, S123) = 12^^ ^[5Ld35.(Si2, S23,pI, S123) + Lc25(Pi, -523, S123)] 

■512-523 

4 I 2 

H ^Lc3(s23, P?, S123) ^Lc2(pi, S23, S123) - 7-^LCi™(Si2, S123) (93) 

S12S23 S12 2S23 

Ld322(si2, S23, pI, S123) = - ( 12 ^^^2 2 ^^ [5Ld35(si2, S23, P?, S123) + Lc25(Pi, -523, Sl23 

^ '5l2'523 

Si2Si23+P?Si23 -2pfsi2 2 x ^12 2m/2 N 

+ S123 2 LC3(S23,Pl, S123) - 7^ LC2 (Pi, S12) 

S12S23 ^-523 

, S12-S123t^ 2m/ \ , ^123 " 2Si2 , 2 \ 

+ S12 ^Tl ^^1 (-^12, Sm) + S123 l^C2(Pi, S23, -5123) 

2S23 -512523 



P1S12 + Sl2£23 — S12S123t „2m^^2 

2si 



+ ^^^""T;^"^"^"^ Lcnp^, S12) ) (94) 



Ld344(si2, S23, pI, S123) = - ( 12 ^^^ 2 ^ [5Ld35(si2, S23, J5?, S123) + Lc25(j5L ^23, S123)] 



'5l2'523 



2P1 ~ 2Si2 ,2 X , 4^1^-512^ / 2 \ 

+ Pl Lc2(j5i,S23,Si23) + Pi 2~LC3(S23, Pi, Sl23j 

S12S23 S12S23 



->2 



+ Si2^i^^Lc^(si2, S123) + ^Lc^ (si2, S123) ) (95) 

ZS23 ^S23 
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B.4 One-mass box 

For box integrals with one massive external leg we use the abbreviation 

Al"' = 2suS23{si23-Sl2-S23). (96) 

The functions we are using for tensor integrals are: 

Ld}™(si2, S23, Su3) = -^^Ldl"^{su, ^23, S123) (97) 



^4 
S12S23 r T ilm/ N , 1 "^123 , 1 "S123- 



Ld2"(si2, S23, S123) = - oAim [si2S23Ld}™(si2, S23, S123) + ^23 In h Si2 In ] (9^ 

<J^4 -523 S12 



Ldr(si2, 523, S123) = -LdJ^(si2, S23, 5123) (99) 

Ld^™(si2, S23, S123) = -— [6Ld^^(si2, S23, S123) + S23Lc?™(s23, S123)] (100) 

S12 

Ld2™(si2, S23, S123) = LdJ'"(si2, S23, S123) - Ld2™(si2, S23, sm)- (101) 
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